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ABSTRACT 

We describe observational evidence for a new kind of interacting-binary-star outburst that involves 
both an accretion instability and an increase in thermonuclear shell burning on the surface of an 
accreting white dwarf. We refer to this new type of eruption as a combination nova. In late 2000, the 
prototypical symbiotic star Z Andromedae brightened by roughly two magnitudes in the optical. We 
observed the outburst in the radio with the VLA and MERLIN, in the optical both photometrically 
and spectroscopically, in the far ultraviolet with FUSE, and in the X-rays with both Chandra and 
XMM. The two-year-long event had three distinct stages. During the first stage, the optical rise 
closely resembled an earlier, small outburst that was caused by an accretion-disk instability. In the 
second stage, the hot component ejected an optically thick shell of material. In the third stage, the shell 
cleared to reveal a white dwarf whose luminosity remained on the order of 10 4 L© for approximately 
one year. The eruption was thus too energetic to have been powered by accretion alone. We propose 
that the initial burst of accretion was large enough to trigger enhanced nuclear burning on the surface 
of the white dwarf and the ejection of an optically thick shell of material. This outburst therefore 
combined elements of both a dwarf nova and a classical nova. Our results have implications for the 
long-standing problem of producing shell flashes with short recurrence times on low-mass white dwarfs 
in symbiotic stars. 

Subject headings: binaries: symbiotic — novae, cataclysmic variables — stars: dwarf novae — stars: indi- 
vidual (Z Andromedae) — stars: winds, outflows — X-rays: binaries 



1. INTRODUCTION 

1.1. Symbiotic-Star Outbursts 

Symbiotic stars are interacting binaries in which ma- 
terial is transferred from an evolved red-giant star to a 
mor e compact, hot s tar, usually a white dwarf (WD: see, 
e.g., iKenvonl 119861: iCorradi. Mikolaiewska. &: Mahonevl 
120031 and references therein). In most symbiotics, the 
red giant under-fills its Roche lobe, and the mass trans- 
fer proceeds via gravitational capture of the red giant's 
wind. An accretion disk may or may not form (depend- 
ing primarily on the relative veloci ty of the re d-giant 
wind and the accreting white dwarf; ILividll988j) . Radi- 
ation from the accreting WD partially ionizes the nebula 
formed by the red-giant wind, and this ionized nebula 
gives rise to optical through far ultraviolet (FUV) emis- 

Smithsonian Astrophysical Observatory, 60 Garden Street, 
Cambridge, MA 02138 

2 NSF Astronomy & Astrophysics Fellow 

3 School of Physics, Trinity College Dublin, Dublin 2, Ireland 

4 School of Cosmic Physics, Dublin Institute for Advanced Stud- 
ies, 5 Merrion Square, Dublin 2, Ireland 

5 Space Telescope Science Institute, 3700 San Martin Dr., Bal- 
timore, MD 21218 

6 The Johns Hopkins University, Department of Physics and As- 
tronomy, 3400 N. Charles St., Baltimore, MD 21218 

7 NOAO, P. O. Box 26732, Tucson, AZ 85726 

8 Astronomy Department, 601 Campbell Hall, University of Cal- 
ifornia, Berkeley, CA 94720 

9 Mullard Space Science Laboratory, University College London, 
Dorking, Surry RH5 6NT 

10 University of Southampton, Southampton S017 1BJ, U.K. 
South African Astronomical Observatory, P. O. Box 9, Ob- 
servatory, 7935, South Africa 

12 NRAO, P. O. Box 0, 1003 Lopezvill Rd, Socorro, NM 87801 

13 UCO/Lick Observatory, Mt. Hamilton, CA 95140 



sion l ines llKenvonllT986l iMu rset et al.lll9911IEsnev 
119961 lEspev fc McCandlissll2003t lNussbaumerll200 

The hot components in symbiotic st ars typically 
have luminosities of around 10 3 L© l)Miirset et alJ 
Il991|) . To produce such a high luminosity by ac- 
cretion onto the WD alone, the accretion rate would 
need to be - 1O -6 M yr" 1 . But it is diffi- 
cult for the WD to accrete from the red-giant wind 
at this rate since red giants in s ymbiotics typically 
lose mass at ~ 10 -7 M© yr" 1 (Scaauist & Taylor 
ll990t). In addition, 10" 6 M© yr" 1 is well above the 
accretion rate needed to produce quasi-steady ther- 
monuclear shell burning iPgc^nsk^^^y^j^^^jS 



Sioji 1 _AcieTnn J _fc TQmczy3 [l£Z9j [Paczvhski ii & ji j^u^als 
"1980 1 jibed 119821 iSion fc Starrfieldl 119861 iSion fc Ready! 
19921 ISion fc Starrfieldl I1994IK ISokoloski et alJ ^200lF 

found that stochastic optical flickering from an accre- 
tion disk, like that seen in cataclysmic variable stars, is 
undetectable in the vast majority symbiotics - probably 
because it is hidden by nuclear-shell-burning light that 
has been reprocessed into the optical by the nebula. In 
fact, the inferred luminosities of ~ 10 3 L© are easily pro- 
duced by the nuclear burning of hydrogen-rich material 
at the rate of a few times 10 -8 M© yr -1 . Accretion onto 
the WD in a typical symbiotic can supply fuel at this rate 
I)liiimall2002j) . Thus,' some level of quasi-steady nuclear 
shell burning on the surface of the WD appears to be a 
common feature of symbiotic systems. 

Symbiotic stars also experience several different types 
of outbursts. The most common type of outburst, 
termed a classical symbiotic outburst, recurs roughly 
every decade or so a nd has optical amplitudes of several 
mag ijKenvod Il986j) . The nature of these eruptions 
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is not yet known. Their peak luminosities appear 
to be too high for t hem to be d i sk instabilities like 
those of dwarf novae l)Kenvon|ll986l iMikolaiewska et al.l 
11995ft . They also recur too frequently to be nova-like 

therm onuclear ru n aways tho se in RS Oph and T 

CrB lEEiS 11982 iFuiimotol 1798% IKenvon fc Trura,nl 
19831 IStarrfield. Sparks, fc Truranl 1985; 

Livio. Prialnik. ^T^^TOvlT989T]Mnk^laiewska et alJl995l) 
or ev en weak shell flashes llSion. Acierno. fe Tomczvkl 
1979). The fundamental difficulty is that most symbi- 
otics contain low-mass WDs. Short recurrence times can 
only be produced for thermonuclear runaways on WDs 
whose masses are very close to the Chand r asekha r limit. 
For example, ISion. Acierno. fc Tomczvkl l)1979f) found 
that even a 1.2-M WD accreting at 10~ 8 M Q yr^ 1 
required more than 600 yr between shell flashes. Thus, 
some additional physics must be involved in classical 
symbiotic outbursts. Since the mass loss during classical 
symbioti c outbursts can take the form of a collimated 
jet (e . g.. iKellogg. Pedeltv. fc Lvonl 120011: iCrocker et ail 
l2001i iBrocksopp et al.l l2004|) . the nature of these 
outbursts is also linked to the issue of jet formation. 

Classical symbiotic-star optical outbursts are tied to 
changes in the hot component of the binary. They can 
therefore be the optical brightening of either the accre- 
tion d isk, the white dwarf, or both l)Kenvon fc Webbinkl 
1984). An accretion disk brightens when mass transfer 
through the disk increases, either because more material 
is being fed to the disk or because of an accretion-disk in- 
stability. The white dwarf can brighten in the optical due 
to either a decrease in effective temperature at constant 
bolometric luminosity, which shifts the WD spectral en- 
ergy distribution to lower energies, or an increase in the 
nuclear burning rate on the white-dwarf surface. A fac- 
tor that complicates this picture slightly is the nebula, 
which reprocesses high-energy photons into the optical, 
where most outbursts are primarily observed. Symbi- 
otic stars can therefore also brighten in the optical due 
to an increase in reproces sed high-energy photons (see 
iNnssbanmer fc Voge]lll98"9l for a model of the reprocess- 
ing). 

Existing models for classical symbiotic outbursts in- 
clude expansion of the WD photosphere at constant 
bolometric luminosity due to the accretion rate ris- 
ing above the maximum value for steady burning 
ijTutukov fc Yun gelsonll9'76llIben lll982l). a thermal pulse 
or shell flash |Kenvcm^TTr>uTanl I198 3 J. or a dwarf- 
nova-like accretion-disk instability i|Duschll Il986albl : 
Mikol aiewska et aTll2002j) . The photospheric-expansion 
model has shown promise for supersoft X-ray sources, 
where it explains the inverse relationship between optical 
and X-ray fluxes seen in some objects (e.g., RX J0513-69; 
iSouthwell et alJll996|) . which arises from changes in the 
location of the peak of the black-body spectral energy 
distribution. Both the thermal-pulse and accretion-disk- 
instability models are related to phenomena that occur in 
cataclysmic variables (CVs) — runaway nuclear burning 
in classical or recurrent novae, and accretion-disk limit 
cycles in dwarf novae. 

1.2. Z Andromedae 

Z Andromeda e (Z And) is the prototype for the class of 
symbiotic stars llKenvonlll986t IMikolaiewska et alJlT988t 
iCorradi. Mikolai^wsk£Tfcl^alion^vn2^03^K It has an 



orbital period of 759 d llFormiggini fc Leibowitd 119941: 
IMikolaiewska fc Kenvonl |l996) , orbital inclination and 
orientation angles of 47 ± 12° and 72 ± 6° respectively 
JSchmid fc Schildl Il997l). and a WD mass of 0.65 ± 
0.28 M ? ijSchmid fc Schildlll997l for a total system mass 
of between 1.3 and 2.3 M ). Distance estimates for Z 
And range from 0.6 kpc t o 2.19 kpc, w ith an average 
value of 1.2 ± 0.5 kpc (see lKennvlll995l and references 
therein) . 

Zandstra-method estimate s of the quiescent W D lumi- 
nosity give Lhot ~ 10 3 L© l)Murset et al.l 11991^ . where 
Lhot is the WD (hot component) luminosity. It is there- 
fore very likely that some degree of nuclear shell burn- 
ing is taking place on the WD surface. Z And is also 
the only symbiotic star for which coherent optical os- 
cillations at the WD spin period of 28 minutes con- 
vincin gly indicate the presence of a strongly magnetic 
WD ( Sokoloski & Bil dstenlll999T) . It was the first binary 
known to contain a WD with both a strong magnetic 
field and quasi-steady nuclear shell burning. Finally, Z 
And is one of a small (^10, but growing) number of sym- 
biotic systems known to produc e collimated outflows i n 
the form of non-relativistic jets ijBrocksopp et al.l f2004'). 

The long-term light curve of Z And sho ws many 
exam ples of classical symbiotic outbursts (see IKenvon! 
1986, p. 4) and demonstrates the non-uniformity of this 
poorly defined type of eruption. The outbursts some- 
times come in pairs (as in 1984 and 1986), or in a se- 
ries of eruptions with decreasing maximum brightnesses 
separ ated in time by slightly less than the orbital pe- 
riod ijKenvon fc Webbinklll984t FFormigg ini fc Leibowitd 
Il994[) . In addition to having different sizes, the outbusts 
in Z And also have different shapes, with some having 
a steeper rise than decline (as in 1997 and comparable 
to dwarf-nova light curves), others having a rise that oc- 
curs at roughly the same pace as the decay (such as in 
1984), and still others having more complex shapes. In 
the 2000-2002 event, the size of the outburst (AV « 2 
mag) was typical for the system, too large to be a simple 
disk instability, and much smaller than outbursts of the 
symbiotic recurrent novae. It was instead closer to the 
~3 mag of a typical symbiotic ( slow) nova, as in AG Peg 
iKenvon. Proga. fc KevesjEOQllh 

iFernandez-Castro et alJ ||1J)95) demonstrated the value 
of multi-wavelength observations of classical symbiotic 
outbursts with their study of the 1984-1986 pair of out- 
bursts of Z And. They found that a shell of material 
was ejected during each of these two outbursts, and that 
the overall behavior of the UV lines and radio fluxes was 
consistent with an accretion-rate fluctuation pushing the 
system over the maximum accretion rat e for qu asi-steady 
nuclear burning. iSokoloskTfc Bildsten (1999), however, 
detected a 28- minute oscillation from magnetic hot spots 
on the WD surface throughout the small outburst in 1997 
and concluded that the WD surface was not hidden by a 
shell ejection during that event. They suggested that it 
could instead have been due to an accretion-disk insta- 
bility. 

To further investigate the nature and cause of classical 
symbiotic-star outbursts, we obtained multi-wavelength 
data, including observations with the FUSE, XMM, 
and Chandra satellites, the VLA and MERLIN radio 
interferometers, and ground-based optical spectroscopy 
and photometry, during the recent 2000-2003 activity 
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Fig. 1. — Long-term V-band light curve from the KAIT monitoring described in t'2. II ('crosses) and the American Association of Variable 
Star Observers (AAVSO; dots), with the times of other observations marked. Ol— 06 refer to selected optical spectroscopic observations 
listed in Table 121 Fl— F9 refer to FUV observations listed in Table El R1-R15 refer to the radio obser vatio ns listed in Tabled X1-X3 refer 
to X-ray observations listed in Table [S] Tl refers to the fast photometric timing search described in H2.5I 



phase of Z And. We describe these observations and 
results in In 33 we use the data to estimate the 
effective-temperature and luminosity evolution of the 
WD, and to investigate the ejection of a shell of material 
associated with the activity in 2000-2001. The WD tem- 
peratures and luminosities help us determine that the 
1997 outburst was due to an accretion-disk instability 
similar to those of dwarf novae, and that the 2000-2002 
event was due to both an accretion-disk instability and 
an increase in the rate of shell burning on the WD. In 21 
we propose a combination nova model for the 2000-2002 
outburst of Z And where the 2000-2003 activity period 
was triggered by an accretion-disk instability, but then 
primarily powered by the resulting increase in nuclear 
burning on the surface of the white dwarf. We summa- 
rize our findings in UT dates are used throughout 
this paper. 



2. OBSERVATIONS AND RESULTS 

Although the hot components (e.g., white-dwarf plus 
accretion disk) in symbiotic stars emit most of their en- 
ergy in the FUV, they also radiate significantly at radio 
through X-ray wavelengths. Important diagnostics are 
found in each of these observational regimes. To make 
progress on the long-standing problem of the fundamen- 
tal nature and cause of classical symbiotic-star outbursts, 
we observed the 2000-2002 outburst of Z And in the ra- 
dio, optical, FUV, and X-rays. A summary of the timing 
and coverage of our data is shown in Fig. d We provide 
a listing of the observations in Tables CHU 



2.1. UBV Nightly Monitoring 

UBV observations were taken nightly starting on JD 
2451802 (2000 September 14) with the 0. 76-m Katzman 
Automatic Imaging Telescope (KAIT; iLi et alJ l2000t 
iFilippenko et al.1 1200H iFilippenkol |2005|) at UCO/Lick 
Observatory on Mt. Hamilton near San Jose, CA. Z And 
is observable with KAIT from about mid-year through 
January, so the KAIT light curves are divided into three 
portions from the latter parts of 2000, 2001, and 2002 
(see Fig. EJl . With integration times ranging from 45 to 
60 s for the U band, 15 to 20 s for the B band, and 5 to 
10 s for the V band, we typically obtained signal-to-noise 
ratios (S/N) of greater than 100 for differential photom- 
etry with respect to two comparison stars in the field of 
Z And. The KAIT light curves plotted in Fig. have 
been converted from the KAIT U BV filter system to the 
standard Johnson system and normaliz ed u sing the abso- 
lute p hotometry of lTomov et alJ l)2003j) and lSkopal et alJ 
( 2002) . We list the results of the photometric monitoring 
in Tabled 

The high quality of the KAIT light curves immediately 
reveals that the rise to optical maximum in 2000 had 
three distinct bursts, separated by two plateau periods. 
The first, second, and third optical rises took approxi- 
mately 2.5, 2.5-3, and slightly more than 3 weeks, re- 
spectively, and the two plateaus lasted for about 1 week 
and then about 1 month respectively. Fig. is a close- 
up view of the first 150 days of the 2000-2002 outburst, 
which clearly shows the three stages. 

Combined with photometry from the American As- 
sociation of Variable Star Observers (AAVSO) and 
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Fig. 2. — Light curves of Z And during the 2000-2002 outburs t. Top panel : full v isual light curve from the AAVSO. Middle panel: 
differential KAIT light curves, normalized to the absolute value s of lTomov et all 1200.^ . The error bars are smaller than the plot symbols. 
Bottom panel: U — B color. Data taken by Sk opal et alj 1200211 before we began our KAIT monitoring are included (large o pen d iamonds) 
to show the blue color spike during the first stage of the optical brightening. The first U — B point, fromlS kopal et all 120021) . is from 
approximately two weeks before the beginning of the outburst. Orbital phase from Mikolaicwska & Kcnvon ( 1996) is plotted on the top 
abscissa. 



iSkopal et al.1 l)2002f) taken just before the beginning of 
the 2000-2002 outburst, the KAIT data show that the 
initial color evolution of the system also had three dis- 
tinct stages (bottom panel of Fig. [21 and Fig. |3J) . Dur- 
ing the first rise, the system became much bluer in both 
U — B and B — V. During the second rise, it became bluer 
still in B — V, but maintained an approximately constant 
U — B color. On the final rise to optical maximum, the 
U — B color again became significantly bluer, suggesting 
a large increase in the temperature and UV luminosity 
of the hot component, whereas the B — V color changed 
very little. 

The optical decline was marked by 0.2-0.3 mag varia- 
tions with no clear pattern. The U — B color remained 
quite blue and fairly constant for roughly the first 300 
days of the decline, and then began to redden. The 
B — V color reddened more gradually throughout the 
entire decline. There was a small rebrightening around 
JD 2452600 (at the end of 2002), where the U - B color 
initially reddened slightly, and then became more blue. 



However, the U — B color never became as blue as when 
the optical flux was at a similar level during the main 
outburst. The peak of the rebrightening occured approx- 
imately 700 days after the main outburst peak. The time 
between the two peaks wa s thus slightly shorter than the 
orbital period of 759 d llFormiggini & Leibowitd [T994t 
Mikol aiewska fc Kenvonll996|) . and similar to that of os- 
cillati ons during the decline of previous large outbursts 
(e.g.. lKenvon fc Webbinldll984T) . 

2.2. Optical Spectroscopy 

P. Berlind, M. Calkins, and several other observers ac- 
quired low-resolution o ptical spectra of Z And with FAST 
ijFabricant et alj[l9 98). a high-throughput slit spectro- 
graph mounted at the Fred L. Whipple Observatory 1.5- 
m telescope on Mount Hopkins, Arizona. They used a 
300 g mm -1 grating blazed at 4750 A, a 3"-wide slit, 
and a thinned Loral 512 x 2688 pixel CCD. These spec- 
tra cover 3800-7500 A at a resolution of 6 A. Starting in 
1994, we obtained over 800 spectra of Z And, with over 
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Fig. 3. — KAIT V-band (crosses), B-band (diamonds), and (7-band (triangles) nightly photometry reveals the stair-step or three-stage 
nature of the rise to optical maximum. V-band data from before the out burst and during the initial fast rise (circled crosses) are from the 
AAVSO. The pre-outburst U and B points are from lSkopal et alj 120021) . 



350 during the 2000-2002 outburst. Several additional 
spectra were obtained by R. Hynes near the optical peak 
with the 4-m William Herschel Telescope (WHT) . We list 
the dates of selected optical spectroscopic observations 
taken near the beginning of the 2000-2002 outburst in 
Table 

We reduced the spectra with IRAF. After trimming 
the CCD frames at each end of the slit, we corrected for 
the bias level, flat-fielded each frame, applied an illumi- 
nation correction, and derived a full wavelength solution 
from calibration lamps acquired immediately after each 
exposure. The wavelength solution for each frame has a 
probable error of ±0.5 A or better. To construct final 
1-D spectra, we extracted object and sky spectra using 
the optimal extraction algorithm APEXTRACT within 
IRAF. Most of the resulting spectra have S/N > 30 per 
pixel. 

Optical spectra of Z And, the prototypical symbiotic 
star, generally have strong high-ionization emission lines 
superimposed on a red-giant absorption spectrum. Our 
spectra show two types of long-term variations. During 
quiescence, high-ionization lines, such as [Fe VII] and the 



Raman-scattered O VI lines at 6830 A and 7088 A, are 
roughly constant with orbital phase. Lower-ionization 
lines, including the H I B aimer lines, decrease signifi- 
cantly in intensity near phase 0.5, when the red giant 
lies in front of the hot component and obscures part of 
the nebula. 

During outburst, all of the lines react to the changing 
luminosity of the hot component as well as the orbital 
geometry. In the 1997 outburst, the high-ionization lines 
increased in intensity relative to the continuum, suggest- 
ing an increase in the effective temperature of the hot 
component. The relative intensities of these lines faded 
during a portion of the 2000 outburst, implying a de- 
cline in the effective temperature of the hot component 
(see H3.1.1I) . Although we do not have high-resolution 
spectra from the 1997 outburst, the WHT spectra taken 
near optical maximum in late 2000 (December 2 and 3) 
show P-Cygni profiles for the He I lines, consistent with 
spectra acquired durin g previous major outbursts (e.g., 
b wings fe Struvdll970T) . At around the same time, FUV 
spectra also show such profiles (see H2.3|) . In Fig. 01 we 
show two example spectra — one from quiescence, and 
one near the peak of the 2000 outburst — both taken just 
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after orbital phase zero. We list the equivalent widths of 

12 lines in 826 spectra taken between 1994 September 

13 and 2003 December 29 in Table El which is available 
electronically. 

2.3. Far Ultraviolet 

Nine observations of Z And (plus one short "safety 
snap") were performed between 2000 November 16 and 
2003 August 4 with the Far Ultraviolet Spectroscopic Ex- 
plorer {FUSE) satellite; see Table H The FUSE spec- 
trographic instrument and on -orbi t performance are de - 
scribed bv lMoos et all l)2000j) and lSahnow et alJ {2000). 
The spectral format consists of eight channels fed by 
four co-aligned prime-focus telescopes and Rowland- 
circle grating combinations, with a separate slit for each 
telescope. In four of the eight channels, LiF overcoated 
optics are used to cover the bandpass from ~ 1000 to 
1190 A. The other four channels use SiC overcoated 
optics to cover the bandpass from the Lyman edge to 
~ 1100 A. Each of the nine observations, consisting of 
anywhere from one to eight orbits, was taken through 
the low wavelength resolution slit (LWRS), which is 
30 x30 square arcseconds in angular area. 

We reduced the FUSE data with CALFUSE pipeline 
versions 2.0.5 and higher, and combined the data for each 
orbit without considering shifts in the stellar position 
within the slits during the Earth-occultation period be- 
tween orbital acquisitions. The individual spectral chan- 
nels were linearized to a 0.01 A scale. They were then 
joined into one spectrum by cross correlation of overlap- 
ping spectral sections. 

For each of the wavelengths listed in TableEJ we calcu- 
lated the continuum level from the unweighted mean of 
all data points at the specified wavelengths ±0.2 A, from 
all the individual exposures at a particular epoch. Wc 
took as the uncertainty the larger of the formal statistical 
error in that mean and the scatter of the measurements 
from the individual exposures. In a few cases, where 
one or two of the individual exposures at a given epoch 
were discordant (centering in the aperture may have been 
poor for several orbits in those cases) , we removed those 
exposures. 

Fig- El shows the joined FUSE spectra (except for the 
safety snap). When the object may have shifted out of 
the slit, we adjusted the flux level to reflect the reduced 
integration time. The spectral continua are rather flat, 
and strong absorption by interstellar H2 spans the entire 
wavelength region (see McCandliss 2003, for H2 absorp- 
tion templates). During the main optical outburst, the 
FUV continuum rose and then fell, peaking at around 
the same time as the optical maximum. The contin- 
uum flux variations are also evident from the 959.5 A, 
1058.7 A, and 1103.4 A flux densities listed in the last 
three columns of Table El By the seventh FU SE ob- 
servation (F7), on 2002 July 5, the FUV flux had re- 
turned to approximately the quiescent value recorded 
by Orpheus and HUT in 1993 and 1995, respectively 
{Orpheus and HUT data were obtained fro m the Multi- 
mission Archive at STScI database, plus see lBirriel et al.1 
119981 ISchmid et al.lll999l) . During the subsequent brief 
optical rebrightening in 2002, the FUV flux dropped by 
almost an order of magnitude, and then rose back to the 
quiescent level again in the final FU SE observation on 
2003 August 4. 



In addition to absorption by some high-ionization 
species, the first FU SE observation shows a large 
amount of cool gas shrouding the FUV continuum. We 
see singly and doubly ionized C, Fe, and Si, as well as ab- 
sorption by P V and the non-resonance 1176 A transition 
in C III. Low-ionization species with P-Cygni profiles in- 
clude N II, P II, C II, and Fe III. The presence of P-Cygni 
profiles in FUSE data is consistent with the optical P- 
Cygni profiles at around the same epoch. By the time 
of the second FU SE observation, the absorbing material 
begins to clear away. Moreover, whereas the first spec- 
trum is almost entirely in absorption, the later spectra 
contain increasingly strong emission features. Toward 
the end of the series of FUV spectra, O VI becomes ex- 
tremely strong as the continuum level falls. 

Fig. shows two examples of the early evolution of 
some line profiles. P V begins in absorption, with a 
blueshifted (P-Cygni-type) line width of approximately 
300 km s . In the second FUSE observation, taken 
immediately before optical maximum, the width of the 
main blueshifted P V absorption decreases to less than 
200 km s _1 , and by the third FUSE observation, taken 
just after optical maximum, the P V absorption compo- 
nents have narrowed further and an emission component 
is present. During this same time, C III A1176 moves 
from absorption to emission. 

The P V absorption-line ratio also demonstrates the 
trend of decreasing opacity evident in the line profiles 
from the first three FUSE spectra. For optically thin 
gas, the ratio of the A1117 to A1128 components of P V 
is 2:1, whereas for optically thick gas, the ratio is 1:1. 
In the first two FUSE observations, the ratio of the 
two P V line components is close to 1:1. By the third 
observation, however, the ratio is approximately 1.33:1, 
indicating that the opacity had begun to decrease. Fur- 
thermore, although the early P V line profiles suggest 
high optical depth, the troughs of the lines do not touch 
zero. Either the r = 1 surface probed in the absorption 
line has significant emissivity relative to the underlying 
continuum, another continuum source is present, or the 
optically thick material only partially covers the contin- 
uum source 14 . The early C III and Si III lines also appear 
to be optically thick, but do not reach zero. 

2.4. Radio 

Radio observations were performed with the NRAO 
Very Large Array (VLA) in New Mexico, and the Multi- 
Element Radio Linked Interferometer Network (MER- 
LIN) in the U.K. Reduction was done with AIPS, and de- 
tails of the reduction and a nalysis of the first 12 obs erva- 
tions (R1-R12) are given in iBrocksoDP et all l|20041. For 
the final three observations (R13-R15; from the VLA), 
we observed in standard continuum mode, with two 
intermediate-frequency (IF) pairs of sidebands record- 
ing both right and left circular polarization, for a total 
bandwidth of 200 MHz at each frequency (100 MHz per 

14 Nebular recombination radiation (predominantly from He II) 
is unlikely to be the possible secondary source of FUV emis- 
sion since the unabsorbed flux expected from the nebula at 
FUV wavelengths is only a few times 10 -13 ergs cm -2 s" 1 
IFernan dcz-Castr o et all 11983. I199fi^ . After extinction by the 
interstellar medium, the nebula contributes on the order of 
10 -14 ergs cm -2 s — 1 at these wavelengths, which is much less than 
the flux we see below the optically thick line profiles. 
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polarization for two polarizations). 

For phase calibration, we used a 110 s / 50 s cycle 
at 8.5 GHz and a 100 s / 50 s cycle at 5 GHz. The 
phase calibrator, 2322+509, is 2.8° away from Z And. 
Flux densities are defined with respect to 0542+498 (3C 
147) and 0137+331 (3C 48), using the 1999.2 VLA flux- 
density scale. The flux densities for Z And represent the 
average of the two IF sideband pairs at each band, in- 
tegrated over a small region surrounding the peak. We 
constructed images of the full primary beam at each fre- 
quency, using natural weighting for maximum sensitivity. 
These images were deconvolved with the CLEAN algo- 
rithm to remove the sidelobes of both Z And and (at 4.9 
GHz) a nearby, extended confusing source. We list radio 
flux densities between 2000 and 2003 in Table EJ The 5 
GHz flux densities are also plotted in the bottom panel 
of Fig. m 

Early in the optical outburst, the 5 GHz flux den- 
sity was very low compared to typical quiescent levels 
of about 1 mjy Our first radio data were obtained ap- 
proximately 40 days after the beginning of the optical 



rise, so the radio flux density could have dropped before 
the initial optical rise, simultaneous with the rise, or af- 
ter the rise. As the optical outburst progressed, Z And 
brightened in the radio and, at 15 GHz, reached 3.5 + 0.1 
mjy, the highest 15 GHz flux density ever recorded for 
this source. The radio brightness peaked between two 
and five months after the optical brightness maximum 
(first at 5 GHz, then at 15 GHz, and finally at 1.4 GHz). 
During the optical rebrightening in late 2002, the radio 
brightness temporarily dropped below typical quiescent 
levels. 

Roughly one year after the start of the optical outburst, 
and about eight and six months after the respective 5 
GHz and 15 GHz brightness maxima, IBrocksopp et alJ 
( 2004) discovered a radio jet from Z And. The extended 
radio component was just 60 mas in length, but was 
oriented in a direction consistent with the outflow be- 
ing almost exactly perpendicular to the binary orbital 
plane Csee lSchmid fc Schildlfl997^ . In 2001 September, 
the ejected component was brighter than the central com- 
ponent. However, the extended component was not de- 
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tcctcd in three subsequent MERLIN observations taken 
over the course of the next eight months. Z And was 
unresolved by the VLA (in C configuration) in 2002 De- 
cember and 2003 January. 

In 2003 July our VLA observations (in A configura- 
tion) indicated that the radio jet was again detectable, 
and had reached ~0.5" to ~ 1" in length. A Gaussian fit 
to the 5 GHz data suggests that the source was extended, 



with a nominal major-axis full-width at half maximum 
(FWHM) of 0.32" ± 0.08" and a minor axis less than 
0.15", at an orientation of 130° ± 15° east of north. 
For comparison, the beam is a Gaussian with FWHM 
0.385" x 0.333" oriented 1.6° degrees east of north. In 
the map, the source shows an elongation to the NW (see 
Fig- 13 j an d we obtained similar results with different 
uii- weighting schemes. The peak flux density of the ex- 
tended component in the 5 GHz VLA image is about 
0.18 mjy beam^ 1 . There is no sign of extended emis- 
sion in the 15 GHz image. The 2003 July VLA extension 
is about 20° west of the direction of the 2001 Septem- 
ber MERLIN jet, but in a similar directio n as a b l ob in a 
1991-1992 VLA+MERLIN map shown bv lKennvl (|1995|) . 

If we take the size of the VLA extended region at face 
value and assume that the material was ejected at the 
start of the optical outburst, the outflow has a velocity 
of around 800 to 1000 km s _1 . The inferred velocity of 
the small (60 mas) jet found with MERLIN was 400 km 
s _1 . However, the two velocity estimates have sufficiently 
large uncertainties that they are roughly consistent. 

2.5. Fast Optical Photometry 

On 2000 November 23, we performed fast optical pho- 
tometry using the 1-m Nickel telescope at UCO/Lick Ob- 
servatory. The goal of these observations was to deter- 
mine whether the 28-minute oscillation due to accretion 
onto the spinning, magnetic WD was detectable during 
the 2 000-2002 outburst, as it w as during the event in 
1997 l|Sokoloski fe Bildsteri ri999L We took repeated 15- 
s exposures using a B filter and an unthinned, phosphor- 
coated, 2048 x 2048 pixel Loral CCD. A typical dead- 
time (with 4x4 binning and fast readout) of 23 s gave 
a time resolution of 38 s. The full time series spanned 
5.6 hours. A description of our differential photometry 
for this field and the data reduction and an a lysis meth- 
ods can be found in ISokoloski fc Bildstenl l|1999li and 
Soko loski et al.1 (|2001^ , respectively. 

The 28-minute oscillation was not detected. After re- 
moval of some hour-time-scale variability with low-order 
polynomial fitting, we found an upper limit on the am- 
plitud e of any 28-minute signal in the data of 3 mmag 
(see lSokoloski et alJ200lL for a description of how we ob- 
tain upper limits using Monte Carlo simulations). The 
conservative upper limit of 3 mmag was obtained with- 
out use of the brightest comparison star, GSC 03645- 
0159 2, since this star is an intrinsically variable ^-Scuti 
star llSokoloski et alJl200l . If GSC 03645-01592 is used, 
and its intrinsic variability removed with a high-order 
polynomial, the upper limit for the 28-minute oscillation 
becomes 1.8 mmag. 

This non-detection indicates that the evolution of the 
28-minute oscillation during the large 2000-2002 out- 
burst was fundamentally different from its evolution dur- 
ing the smaller 1997 outburst. In 1997, the oscilla- 
tion amplitude changed roughly in step with the hot- 
component optical flux, reaching almost 6 mmag when 
the optical flux was 0.9 mag above the pre-outburst level. 
In 2000, despite the fact that Z And had brightened by 
1.8 mag compared to quiescence, the oscillation was not 
detected. 

2.6. X-Ray 




Fig. 7. — Extended 5 GHz radio emission detected on 2003 July 24 with the VLA (in the A configuration). The extension is in a direction 
roughl y 20° west of th e small, transient ejection seen in 2001 September, just one year after the beginning of the optical outburst, by 
IBrocksoPD et alj <200ll) . The contour levels are 6.1 X 10~ 2 X (-4, -2.83, -2, 2, 2.83, 4, 5.657, 8) mjy. 



We obtained three X-ray observations of Z And during 
the first year of the 2000-2002 outburst. The first ob- 
servation (XI ) was performed with the Chandra X-Ray 
Observatory {Wcisskopf et aD l2002j) on 2000 November 
13, just over two months into the optical outburst, at the 
beginning of the third stage of the optical rise. The sub- 
sequent two X-ray observations (X2 and X 3) were per- 
forme d with the XMM-Newton Observatory ijAschenbachl 
12002ft . on 2001 January 28 (~50 d after the optical out- 
burst peak) and 2001 June 11 (6 months after the opti- 
cal peak), respectively. See Table El for a summary of all 
three X-ray observations. 

Observation XI was performed with the Chandra High 
Energy Transmission Grating Spectrometer (HETGS). 
The data were processed using updated calibration files 
(CALDB 2.27), and reduced with the software package 
CIAO v3.1. In 19 ksec, the source was detected with 109 
source counts in the zeroth-order spectrum. Due to the 
low number of counts, we used only the zeroth-order data 
for our analysis. To reduce the instrumental background, 
we only selected photons with energies between 0.3 and 7 
keV. Based on examination of source-free regions in the 
SI chip, no periods of high background occurred during 
the observation. 



We performed spectral analysis for all X-ray data using 
the software packages Sherpa (within CIAO) and XSPEC 
vll.3. For observation XI, we extracted the source spec- 
trum from a 3"-radius circular region centered on the 
source, and used an annular region also centered on the 
source for the background. We rebinned the spectrum 
to have at least 10 counts per spectral bin. In the en- 
ergy range over which the source was detected, the back- 
ground contributed less than 0.1%. 

Observations X2 and X3 were performed with the 
XMM-Newton European Photon Imaging Cameras 
(EPIC), in full- frame mode, using the medium neutral- 
density filter. After rejecting data with high background 
levels, we retained good-time intervals of 16 ksec and 
4 ksec (out of 20 ksec and 15 ksec) for observations X2 
and X3, respectively. 

We reduced and analyzed the XMM-Newton data with 
the software package SAS v5.4.1, using only data in the 
0.2-10 keV range. We extracted the EPIC-pn source 
spectra from a 30"-radius circle centered on the source, 
and the background from a source- free region. Z And was 
detected in the EPIC-pn with over 1500 source counts 
during observation X2, but was significantly fainter 
(about 140 EPIC-pn source counts) during observation 
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Fig. 8. — Left: X-ray spectrum of Z And from the 2000 November Chandra observation (XI), fitted with an absorbed blackbody model 
(kT = 0.2 keV) plus an absorption edge at ~ 1 keV. Right: X-ray spectrum of Z And from the 2001 January XMM-Newton observation 
(X2), fitted with an absorbed blackbody (kT = 0.11 keV) and powerlaw (a = 0.79) model with absorption edges at 0.29 keV (most likely 
due to a calibration error), 0.64 keV, and 0.96 keV. 



X3. In addition, the background level was significantly 
higher during observation X3 (approximately 0.7 counts 
s _1 , probably due to solar activity, compared to < 0.05 
counts s' 1 in observation X2). To improve the statistics, 
wc rcbinncd the spectra from X2 and X3 to have 20 and 
10 counts per spectral bin, respectively. 

The X-ray emission detected in observation XI fell pri- 
marily between 0.3 and 2 keV. With only ~100 source 
counts and a large uncertainty in the amount of intrin- 
sic absorption, we are unable to use the spectrum from 
XI to strongly constrain physical models. Neverthe- 
less, we fit the spectrum with several different models 
as a means of parametrizing its softness and investigat- 
ing the absorption. We could not obtain an acceptable 
fit for a single-component Raymond-Smith (RS) plasma 
model. For a powerlaw model, the best fit is statisti- 
cally acceptable, but gives an unphysically large pho- 
ton index of T = 5.9 (where the photon flux is given by 
A(E) = KE~ T , with E being energy and K a normal- 
ization constant). The absorbing column density in the 
powerlaw model is Nh = 5 x 10 21 cm~ 2 , which is sig- 
nificantly hig her than the Gala ctic value of 1.4 x 10 21 
cm~ 2 (Di ckev fc Lockmanl lT990). but lower than values 
of ~ 10 22 cm~ 2 sometimes found for symbiotic stars near 
orbital phase zero (where X-ray emission from the accret- 
ing white dwarf must pass throug h the maximum amount 
of the red-giant wind: lVogell99lD . As discussed in M3.1.2I 
we adopt E(B — V) = 0.27 mag for Z And. Based on the 
re lationsh i p betw een E{B — V) and Nh shown in Fig. 1.1 
of lSpitzerll|1978fl . E(B~V) = 0.27 m ag is consistent with 
the Ga lactic column density given bv lDickev fc Lockmanl 
lH990l) . 

A blackbody model yields a good fit if we include an 
absorption edge at ~ 1 keV; the significance level of the 
additional absorption edge is greater than 95%. The 
resultant blackbody temperature is 0.21±0.04 keV, or 
T bb = (2.4 ± 0.5) x 10 6 K, and in this model N H is con- 
sistent with the Galactic value. The absorbed flux in the 
energy bands 0.3-7 keV and 2-20 keV is (2.5±0.2) x 10~ 13 
erg cm -2 s _1 and (6.9 ± 0.7) x 10~ 15 erg cm~ 2 s -1 , re- 
spectively. The spectrum from observation XI, with the 
blackbody-plus-absorption-edge model, is shown in the 



left-hand panel of Fig. [S] We list the best-fitting spec- 
tral parameters in Table 

For observation X2, there is no single-component spec- 
tral model which can provide an acceptable fit to the 
data, in part because low-level emission was detected 
all the way out to 10 keV. As with observation XI, 
a single-component RS plasma model does not gener- 
ate a good fit. Instead, to obtain a reasonable xl> a 
blackbody component with kT = 0.1 l±8l8i keV ( T <* = 
(1.3±0.2)xl0 6 K), a powerlaw component with T = 0.79, 
and two absorption edges (at 0.64 and 0.96 keV) are re- 
quired. The optical depth of both edges is ~1. An ad- 
ditional edge at 0.29 keV is probably due to the calibra- 
tion error near the carbon edge at 0.3 keV. The derived 
Nh in this model is consistent with the Galactic value. 
The powerlaw component comprises about 6% of the to- 
tal unabsorbed flux. The absorbed flux in 0.3-7 keV 
and 2-20 keV is (1.3 ± 0.04) x 10" 13 erg cm" 2 s" 1 and 
(6.8 ±0.2) x 10~ 14 erg cm -2 s , respectively. The spec- 
trum from observation X2 and best-fitting spectral model 
are shown in the right panel of Fig. [S] The spectral-fit 
parameters are listed in Table [7| 

Z And was significantly fainter during observation X3. 
With the low source count rate and high background level 
during this observation, it is clear that the spectrum was 
soft, but it is otherwise not tightly constrained. The 
spectrum is best described by an absorbed blackbody 
model with kT = 0.12 ± 0.07 keV (T bb = (1.4 ±0.8) x 10 6 
K), and N H = 2.4^ x 10 21 cm" 2 A powerlaw with 
an unphysically high photon index of > 10 produced a 
similar-quality fit, as did a RS plasma model fit, which 
gave kT = 0.086 ±0.002 keV (T RS = (1.00 ±0.02) x 10 6 
K) and an absorbing column close to 10 22 cm -2 (Nh = 
(9.5 ± 0.05) x 10 21 cm' 2 ). The absorbed fluxes in 0.3-7 
keV and 2-20 keV are (2.0 ± 0.5) x 10~ 14 erg cm- 2 s' 1 
and (6.1 ± 2) x 10~ 18 erg cm -2 s _1 , respectively. 

In all three X-ray observations, we find formally ac- 
ceptable fits to blackbody models with temperatures of 
a few times 10 6 K. However, these temperatures are 
too high for the X-ray emission to be blackbody emis- 
sion from either the WD surface or accretion disk. A 
WD with a radius of 10 9 cm and blackbody temperature 
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Tbb = 2 x 10 K would have a luminosity of 3 x 10 Lq 
— two orders of magnitude greater than the Eddington 
limit for a 0.65 M© WD. Moreover, our optical and FUV 
data indicate that the effective temperature of the WD 
in Z And typically does not rise above 180,000 K (see 
H3.1(l . Given that we have very poor statistics for two 
of the three X-ray observations, that symbiotics gener- 
ally have rich line spectra at other wavelengths, and that 
there are multiple potential sites of X-ray emission and 
absorption in Z And, we suspect that a more complex X- 
ray spectrum, possibly containing overlapping emission 
lines from a plasma and absorption features from the 
ionized nebula, is mimicking a blackbody in our CCD- 
resolution spectra. 

In both the FUV (see and radio (see fl2~4) data, 
we find evidence that a significant amount of material 
was ejected from the WD into the dense symbiotic neb- 
ula. Thus, an alternative possible source of X-ray emis- 
sion from Z And during outburst is shock- heated plasma. 
Taking a distance of 1.2 kpc, the unabsorbed X-ray fluxes 
from Table [7| give 0.3-7.0-keV X-ray luminosities of al- 
most 10 32 erg s" 1 for observation XI, dropping to a few 
times 10 31 erg s -1 by observation X3. These X-ray lu- 
minosities are more than five orders of magnitude less 
than the total bolometric luminosity of the WD near the 
peak of the outburst (see fl3.1.2|) . Thus, most of the out- 
burst energy is emitted at slightly lower than X-ray ener- 
gies, as expected given the WD effective temperature of 
< 180, 000 K. The total 0.3-7-keV X-ray energy radiated 
in four months (the approximate time for which the X- 
ray emission was high compared to quiescence) is about 
8 x 10 38 erg — less than the kinetic energy of 10 _9 M© of 
material moving at 300 km s _1 (the approximate velocity 
of the outflowing material; see § 12.31 and § 12.4ft . Since at 
least 10 _7 M© was probably accreted in 2000 (see JP, the 
availability of this much kinetic energy for conversion to 
X-ray emission is plausible. In addition, the temperature 
behind a strong shock moving at 400 km s^ 1 is roughly 
2 x 10 6 K (10 6 K if the shock velocity is 300 km s" 1 ), 
and plasma at this temperature will produce soft X-ray 
emission. Moreover, if we take a cooling time of a few 
months and a plasma temperature of 10 6 K, we infer a 
post-shock density on the order of 10 7 cm -3 . This esti- 
mate gives a pre-shock density of a few times 10 6 cm~ 3 , 
which is reasonable for the outer nebula in Z And (as- 
suming that the density d ecreases from the inner-nebula 
value of ~ 10 10 cm~ 3 , from lFernandez-Castro et alll988l 
as the inverse of the distance squared) . 

The presence of absorption edges provides additional 
evidence for significant amounts of hot, possibly shock- 
heated gas. In observation XI, we found a 1.02 keV edge 
that is close to the energies of either the Fe VII Li edge at 
1.03 k eV or the Ne VI K edge at 1.04 keV. lEbisawa et al.l 
l)2001l) saw the same edge in observations of the supersoft 
X-ray source RX J0925. 7-4758, and they ascribed it to 
a high column of ionized material between the WD and 
the observer. We also detected a 0.96 keV edge in the 
observation-X2 spectrum, possibly related to the same 
feature. The 0.64 keV edge seen in the observation-X2 
spectrum is close to the Lyman edge of N VII at 0.67 keV. 
A similar feature was fou nd in the X-ray spec trum of 
the recurrent nova U Sco fKaha bka et"atl fl999V since 
it is hard to imagine that the hot-component effective 
temperature could get high enough to produce N VII via 



photoionization in Z And, the presence of the 0.64 keV 
edge could be evidence for shock heating. The hard tail 
in observation X2 provides additional evidence for shock- 
heated gas at this time. 

Comparing our X-ray fluxes to those from a 

1993 quiescent-sta te ROSAT PSPC observation 
ijMiirset et alJll997[K the X-ray flux of Z And increased 
by about an order of magnitude during the 2000-2002 
outburst. The PSPC count rate in the ROSAT energy 
band (0.12-2.4 keV) in 1993 was 0.003 ± 0.002 counts 
s _1 , equal to approximately 3 x 10~ 14 ergcm~ 2 s _1 
(for reasonable blackbody spectral models). In obser- 
vations XI and X2, the 0.3-2.0 keV flux was at least 
10 -13 ergcm~ 2 s _1 or a few times 10 -13 erg cm~ 2 s _1 . 
By the time of observation X3, six months into the 
optical decline, the X-ray flux had decreased to a level 
similar to that measured by ROSAT during quiescence. 

We also ex amined two public Rossi X-ray Timing Ex- 
plorer PCA l).Iahoda et alJll996T) observations of Z And 
— one during the decline of the small optical outburst 
in 1997 (45 ksec; 1997 November 28), and another a year 
later after the system had returned to quiescence (12 ksec 
and 37 ksec on 1998 November 17 and 18). In the first 
RXTE observation, Z And was detected with a rate of 
~2 counts s _1 . Using the software PIMMS and assuming 
a powerlaw spectral model, this count rate gives a 2-20 
keV flux of a few times 10 -12 erg cm -2 s , two to three 
orders of magnitude higher than the 2-20 keV fluxes es- 
timated for the 2000-2002 outburst. Z And was not de- 
tected during the second, quiescent-state RXTE obser- 
vation. Thus, unlike in the supersoft X-ray sources, some 
observations of which show an inverse relationship be- 
tween the optical and X-ray brightness ((Southwell et alJ 
1996), in Z And the X-ray flux rose during both the 1997 
and 2000-2002 optical outbursts. 

3. ANALYSIS AND CONCLUSIONS 

To investigate the changes in the hot component from 

1994 to 2003, and to diagnose the nature of the two out- 
bursts during this time period, we use the data described 
above to estimate the effective temperature and luminos- 
ity of the WD-plus-disk as a function of time. In addi- 
tion, optical and FUV line profiles provide information 
about material ejected during the 2000-2002 outburst. 

3.1. Hot- Component Effective-Temperature and 
Luminosity Evolution 

3.1.1. Effective Temperature 

We estimated the effective temperature of the hot com- 
ponent between 1994 and 2003 usi ng both the ratio of 
nebular emission lines l)Iiiimalll981[) and the presence of 
different ionization states in th e optical and FUV spectra 
l)Miirset fc Nussbaumerl ll994jL Assuming a radiation- 
bounded nebula (all pho tons capable of ionizing H and 
He ionize these atoms), lliiimal l(198lD showed that the 
effective temperature of a central source of ionizing pho- 
tons is 

^•° i8 - 3 K ii6jgfs ) 1/2+5 - 13 - (i » 

where -F4686 , Fh0, and -F4471 are the nebular emission- 
line fluxes of He II A4686, H/3, and He I A4471, respec- 
tively. This expression is valid for effective temperatures 
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between 70,000 K and 200,000 K. In our treatment, we 
neglected the A4471 flux, because F4471 < 0.1 Fh/3- To 
minimize calibration errors, we used the ratio of equiv- 
alent widths, EW4686/EWff,3, instead of the flux ratio. 
Because the continuum is nearly flat at AA4650-4900, 
this approximation is justified. Errors in the equivalent 
widths are ~ 5%-10%; thus, the statistical error in the 
effective temperature is ~ 10%— 15%. To correct for the 
fact that near orbital phase zero the H/3-emitting region 
is partially occulted by the red giant, and thus the H/3 
equivalent widths (EW^) do not reflect the true H/3 flux 
from the nebula at these times, we multiplied Tnji ma by a 
sinusoidal factor, Thot = Tuji ma (4/ (5+cos(2ir0))), where 
9 is the orbital phase. This expression was empirically 
derived to force the quiescent orbital-phase-zero values of 
Thot on average to equal the quiescent orbital-phase-0.5 
values. 

We also estimated Thot based on the presence of dif- 
ferent ionization states in the optical and FUV spectra. 
iMiirset fc Nussbaumerl l)1994t) derive a relation between 
Thot (in degrees K) and Xmax, the ionization potential 
(in units of eV) of the highest-ionization species visible 
in the spectrum: 

T ho t = 1000x max ■ (2) 

They used photoionization models to estimate an accu- 
racy of ~10% for this expression. The ionization-state- 
based estimate of Thot for Z And generally agreed with 
the line-ratio-based estimate described in the previous 
paragraph. However, the two methods gave slightly dif- 
ferent results at the time of maximum WD luminosity 
during the 2000-2002 outburst, when the assumption 
of a radiation-bounded nebula might have been violated 
(the rather low X-ray absorpting column and the higher- 
than-average radio flux are both consistent with the neb- 
ula being more ionized around the time of peak lumi- 
nosity), and near orbital phase zero in 2002, when the 
H/3-emitting region was partially occulted (as described 
above). Because of these problems with the line-ratio- 
based estimates, and because the combination of optical 
and FUSE spectra provides a good sampling of ioniza- 
tion states that bracket the temperatures of interest, the 
ionization-state-based results were given precedence on 
these occasions. 

For our spectra, the persistence of He II A4686 (x = 
54.4 eV) indicates T hot > 55,000 K throughout 1994- 
2003. Both Raman-scattered O VI at A6830 (x = 
114 eV) and the [Fe VII] A6087 (x = 100 eV) lines 
are also generally visible when Z And is inactive, sug- 
gesting Thot > 10 5 K in quiescence. During the 1997 
outburst, Raman A6830 strengthened considerably, con- 
firming the increase in 7} lot to 180,000 K suggested by 
F4SSg/Fh/3- During the rise to optical maximum in 2000, 
[Fe VII] A6087 and Raman A6830 weakened and then 
disappeared, and in FUSE spectra, O VI AA 1032,1038 
was absent, all supporting the decline in Thot to ~ 
80,000-90,000 K. These lines reappeared during the op- 
tical decline, confirming the rise in effective temperature 
determined from F^ge/Fiip- In FUSE spectra F4 - 
F6, taken midway through the optical decline in 2001, 
O VI AA 1032,1038 are strong, but Ne VI AA 1006,1010 
(x = 126 eV) are weak or absent. For this time period, 
we adopt 120, 000 K, slightly below our F^ s& /Fhp esti- 
mates. In FUSE spectra F7 and F9, in 2002 July and 



2003 August, Ne VI AA 1006,1010 are present, confirm- 
ing the increase in Thot back to the pre-outburst value, in 
accord with the simultaneous return to the optical quies- 
cent state. During the brief optical rebrightening in late 
2002, the fading of Ne VI AA 1006,1010 on FUSE spec- 
trum F8 indicates that Thot temporarily decreased, and 
we again adopt a value of 120, 000 K for this time period. 
Overall, our line-ratio-method and ionization-state-based 
estimates of Thot appear to be consistent to within 10- 
20%. The hot-component effective temperature, Thot-, is 
plotted as a function of time in the top panel of Fig. 

We found that T ho t = 150, 000 ± 7, 000 K during qui- 
escence. Our value of the quiescent WD effective tem- 
perature is somewhat higher than the quiescent-state 
Za nstra temperature o f 120,000 K (on average) found 
by IMiirset et all (|1991l) between 1979 and 1987 and o f 
roughly 10 5 K found by iFernandez-Castro et all (|198ST> 
between 197 8 and 1982. Howeve r, the effective temper- 
atures from IMiirset et af] l|1991fl increase during 1979- 
1981 (108,000 K on average) and 19 84-1987 (126,000 K 
on average), and lBirriel et all l|1998[) find T = 111,000± 
4, 000 K from a FUV observation in 1995 March. Our 
higher values could thus refl ect long-term chang es of 
a similar magnitude in Thot- iTomov et alJ f:2003) find 
Thot = 150,000 K, as we do, in 1999 September. 

In summary, the effective-temperature history of the 
WD is markedly different for the 1997 and 2000-2002 
outbursts. During the 1997 outburst, Thot increased by 
about 20%, to more than 180,000 K. It returned to the 
quiescent-state value by midway through the optical de- 
cline. Thot showed a more complex evolution durin g the 
2000-2002 event. Photometry from lSkopall l)2003j) indi- 
cates that very early in the 2000-2002 event (within the 
first few days), the [/-band flux jumped by 2 mag and 
the U — B color decreased (i.e., became bluer; see ISkonall 
120031 for more complete coverage of this blue spike), sug- 
gesting that Thot initially increased. Within two weeks of 
the start of the outburst, however, as the optical flux con- 
tinued its rise to maximum, the Thot evolution changed 
direction. Over the next two months, Thot decreased to 
< 90,000 K. By the time of optical maximum, Thot was 
beginning its recovery. It returned to the quiescent-state 
value before the completion of the optical outburst. Fi- 
nally, Thot again dropped briefly during the small optical 
rebrightening in late 2002. The dominant behavior of 
Thot dropping as the optical flux rises to maximum, as in 
the 2000-2002 event, is typical of historical outbursts. 

3.1.2. Luminosity 

To determine the bolometric luminosity of the 
WD throughout the 2000-2002 outburst, we esti- 
mated the WD radius at the time of each FUSE 
ob servation by s c aling the WD photosphere models 
of iBarman et alJ l)2000|) to the extinction-corrected 
FUSE fluxes, using the effective temperatures de- 
scribed ab ove, a distance of 1.2 kpc, and a mass of 
0.65 M Q ISchmid fc Schildl Il997l) . V alues of E(B - 
V) in the literature in clude 0.3 mag jAjtamor e et all 
rmsT^I . 0.27-0.29 mag (IKenvori fc Wehhmld 11 9841) . 0.3 
mag llMikolaiewska fc Kenvonl Il996|) . and 0.35 mag 
l|Viotti et al.l 11 982ft . We take FAB - V) = 0.27 mag, 
the lowest of these values. This choice is the most 
conservative for our purposes, as larger values will pro- 
duce even greater WD bolometric luminosities. Further- 
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Fig. 9. — Top panel: the hot-component effective temperature, T^ot, between 1994 and 2003. The large blue crosses mark times when our 
two methods for estimating T^ ot disagreed slightly, and the ionization-state- based estimate was used (see text for details). T^ ot rose during 
the 1997 outb urst and drop ped during the 2000-2002 outburst. f/-band fluxes and U — B colors near the very beginning of the 2000-2002 
outburst (from Skopal 2003) indicate that Tft oi probably rose briefly before dropping in 2000 (red arrow). Middle panel: V-band light curve, 
from the AAVSO, showing both the 1997 and 2000-2002 outbursts. Bottom panel: 5 GHz radio flux densities, which drop (compared to 
the ty pical quiescent level of about 1 mjy) early in the 2000-2002 outburst. Orbital phase from the ephemeris of lMikolaiewska fe Kenvonl 
1 1996) is shown at the top of the plot. 



more, this low value of E(B — V) produces a dereddened 
flux ratio (-FA1059/-FA1103) that is more consistent with 
the theoretical expectation, and also more consistent 
with the interstellar reddening from dust em ission maps 
l|Schlegel et alJll998t iBurstein fc Heilesll 1982ft . With the 
WD photosphcric radii and effective temperatures deter- 
mined previously, we were then able to obtain the lumi- 
nosities using the standard relation Lhot — ^Rhot^hot- 
The FUSE fluxes provide a good approximation to the 



flux directly from the WD because the contribution from 
both the red giant and the nebula are negligible in the 
FUV (for a complete spectral energy d istribution and 
model, see iFernandez- Castro et aLll!988ft . 

The derived hot-component luminosities, which we list 
in Table |3 are close to 10 4 Lq from the end of 2000 to 
late 2001. Although these results are quite sensitive to 
the reddening used to correct the FUV fluxes for extinc- 
tion, our choice for E(B — V) on the low end of the range 
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of published estimates means that the true luminosities 
could have been even higher. In addition, we find rea- 
sonable agreement between our estimate of Lhot ~ 1500- 
2000 Lq in quiescence (FUSE observations F7 and F10) 
and the quiescent-stat e Lhot estimates of a pproximately 
620-1600 L Q found bv lMiirset et alJ l|199lD . who used a 
modified Zanstra method, which is much less sensitive to 
reddening. For comparison with our outburst luminosi- 
ties, the Eddington luminosity for a 0.65 Mq white dwarf 
is L Ed d = 3 x 10 4 L© (M/0.65M©). 

3.2. Disk Instability in 1997 

We interpret the 1997 outburst in terms of a disk in- 
stability that produces a sudden increase in the accretion 
rate through the disk surrounding the white dwarf. Con- 
siderable evidence supports this interpretation: 1) the 
amplitude and shape of the optical light curve; 2) the 
behavior of the short-period photometric oscillation due 
to magnetic accretion onto a rotating white dwarf; and 
3) the evolution of Thot m 1997. 

Given the small size and short duration of the 1997 
outburst, it was unlikely to have been caused by a 
recurrent-nova-like thermonuclear runaway or a weak 
shell flash. Thus, the main alternative to the disk- 
instability interpretation is an expansion of the white- 
dwarf photosphere at roughly const ant bolometric lumi- 
nosity jT u^ukov Yung elsonlll976T) . In symbiotic stars, 
iKenvon fe Webbinkl l)1984f) showed how the spectroscopic 
evolution of this alternative differs from the evolution ex- 
pected from a disk instability. Here, we consider each of 
the pieces of evidence listed above and use this informa- 
tion to reject the expanding-photosphere picture in favor 
of a disk instability. 

The triangular shape of the light curve during the 
1997 outburst r esembles the outbursts of many dwarf 
novae (see, e.g., lWarnerlll995|) . The rapid rise time of 
< 20 days is similar to the typical rise times of ~ 2-10 
days for dwarf novae, especially given that longer dwarf- 
nova rise times are typically seen in CVs with longer or - 
bital periods (and therefore larger disks; Warner 1995). 
The slower, nearly linear decline is also reminiscent of 
dwarf novae, where a cooling wave propagates inward 
through the disk at the local sound speed. In dwarf no- 
vae, the cooling wave takes, for example, 7-10 days 
to travel roughly 0.7 R©. The longer 300 day decline 
in Z And is reasonable for a disk with an outer radius 
of 20 - 40 R©, or roughly 3 - 7% of the binary sepa- 
ration (10 - 20% of the hot-component Roche-lobe ra- 
dius if we take Mwd — 0.65 Af© and a red-giant mass 
of 1.3 M©). Moreover, the faster rise than decline is 
suggestive of an outburst in which the heating wave is 
initia ted at large radii and moves inward through the 
disk ijCannizzo. Wheeler, fc Polidanl 1 1986(1 . Besides the 
shape, the amplitude of the 1997 outburst is close to the 
~1 mag maximum amplitude expected for a disk insta- 
bility in a symbiotic star (the amplitude is low because 
of the additional light from the red giant and nebula; 
IKeuvonllT98?i IPnschlllT986hl) . 

Dur ing the small 1997 outburst, ISokoloski fc Bildstenl 
( 1999) detected a 28-min photometric oscillation when 
the system was near optical maximum and then through- 
out the decline back to optical quiescence. They inter- 
preted this oscillation as a hot-spot on the surface of 
a white dwarf with a rotational period of 28 min. The 



properties of the light curve are consistent with a hot spot 
produced by accretion channeled from the disk to the 
white dwarf by a magnetic field with a surface strength 
of at least 3 x 10 4 G (or 6 x 10 6 G if the WD is in spin 
equilibrium). The behavior of the oscillation through- 
out the outburst provides a clear discriminant between 
a disk instability and an expanding photosphere. In the 
expanding-photosphere model, the optical emission from 
accretion remains constant while optical emission from 
the WD photosphere increases. Therefore, the ratio of 
optical flux from accretion to that from the WD pho- 
tosphere drops by at least the same factor as the over- 
all increase in optical brightness. Thus, the expanding- 
photosphere model predicts that the pulse fraction of the 
oscillation will decrease during an outburst. In contrast, 
because the luminosity of an accretion hot-spot is likely 
to track the overall accretion luminosity, the pulse frac- 
tion of the oscillation should remain roughly constant 
during a disk instability. 

Our data clearly favor a disk instability as the source 
of the 1997 outburst. The pulse fraction of ~ 0.6% near 
the peak of the outburst is comparable to the pulse frac- 
tion of ~ 0.7% during quiescence (after subtracting out 
the approximate contribution of the red giant to the B- 
band flux, which we take to be 75% in quiescence; see 
Fcriiandcz-Cas tro et alJ Il988[) . In contrast, we would 
expect at least a factor of three reduction in the pulse 
fraction for a 1 mag rise in brightness due to an expand- 
ing photosphere. The data clearly rule out such a large 
change in the pulse fraction. 

The derived change in Thot also favors a disk instability 
over an expanding photosphere. In the Rayleigh-Jeans 
tail of the WD blackbody distribution, the optical inten- 
sity is proportional to Thot- The optical flux is therefore 
proportional to R\ ot Thot- If the bolometric luminosity 
is constant, as is the case in the expanding-photosphere 
model, then B\ ot oc T^ 4 , and the optical flux is pro- 
portional to Tj^.. For a 1 mag optical eruption, the 
expanding-photosphere mod el thus predicts a factor ~ 
1.3 decline in T hot (see also lBath fc Shavivl H976). Our 
analysis of the optical emission lines indicates a factor of 
1.2 increase in Thot (Fig. EJ- The disk-instability model 
is clearly more consistent with the data. 

3.3. Disk Instability Plus Enhanced Shell Burning in 
2000-2002 

The hot-component effective-temperature evolution 
was quite different during in the 2000-2002 outburst. 
These two eruptions were therefore examples of distinct 
types of events. Unlike the 1997 event, the properties of 
the 2000-2002 outburst are not consistent with a simple 
disk-instability model. Evidence against the simple disk- 
instability model for this eruption includes the following: 
1) the rise to optical maximum proceeded in three dis- 
tinct stages; 2) the optical brightness increased by well 
beyond the 1 mag th at can be prod uced by a disk insta- 
bility in a symbiotic l)Kenvonlll986f) : 3) Thot evolved in a 
rather complex way; 4) the optical oscillation at the WD 
spin period due to magnetic accretion was absent near 
the peak of the outburst; and 5) the amount of material 
needed to power this event by accretion alone is much 
higher than could reasonably have been accreted during 
the outburst. 

The start of the 2000 eruption closely resembled the 
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evolution of the 1997 outburst. Fig. llQI in which the 1997 
and 2000 light curves are overlaid, shows the similarity 
between the first stage of the 2000-2002 event and the 
rise to maximum in 1997. Very early in the 2000-2002 
event, Thot appears to have increased, as it did in 1997 
(see j jSZDJ ■ We do not have the spectral data needed to 
determine Thot during the first week or so of the 2000- 
2002 outburst, or the fast photometry to show that the 
pulse fraction of the oscillation remained constant dur- 
ing the early rise in 2000. However, the U — B color 
evolu tion during the first week of the outburst ijSkopall 
2003) strongly suggests that Thot increased during the 
early rise. Furthermore, since the U brightness tends to 
be dominated by reprocessed ionizing radiation f rom the 
hot c omponent, the early jump in [/-band flux ijSkopall 
l2003|) indicates that the 2000-2002 outburst was initiated 
by a process that caused the luminosity of the hot com- 
ponent to increase rapidly. We therefore conclude that 
the 2000-2002 outburst began with a disk instability, as 
in 1997. 

If the 2000-2002 outburst was entirely accretion pow- 
ered, the required accretion rate to produce such a high 
luminosity would be 

L Rhot 



GAL 



WD 



10 4 £ m 0.1 Rr. 



AL 



WD 



nrr W , M G yr" 1 , (3) 
0.65 Mq J 

where L is the bolometric luminosity, Rhot is the WD 
radius, Mw n is the WD mass, and Rhot = 0.1 Rq as 
indicated bv lMiirset et al.l l)199lH and our scaling of the 
FUV fluxes to photospheric models (see Table It is 
difficult to imagine how such a high accretion rate could 
be sustained in Z And for a full year. 

If the outburst was nuclear powered, on the other hand, 
the burning of only a few times 10 -7 Mq of fuel can pro- 
duce Lhot ~ 10 4 Lq for one year. Moreover, since Lhot 
was higher by an order of magnitude in outburst com- 
pared to quiescence, and we argued in § II .11 that the qui- 
escent hot-component luminosity already indicates the 
presence of nuclear shell burning, it follows that the rate 
of nuclear burning must have increased during outburst. 
iKilpio et al.l l)2005|) also suggest that this outburst in- 
volved enhanced nuclear burning (but as a result of the 
collapse of the disk when the wind speed dropped be- 
low a critical level). If we include the kinetic energy of 
ejecte d mass in the total energy produced by th e out- 
burst l|BrocksoDD et all 12004 iTomov et all 120031. ther- 
monuclear involvement is even more strongly indicated. 

3.4. Shell Ejection and Dissipation 

Our observations suggest that during the second stage 
of the 2000-2002 outburst, a shell of material (or pos- 
sibly an optically thick wind) was blown from the sur- 
face of the WD. During the third stage, this ejected 
shell became optically thin. During the second stage, 
Thot dropped (see Fig. |5J) and the U — B color reddened, 
first sharply, and then more ge ntly (see Fig. EJ. The 
FUSE spectrum taken near the end of the second stage 
was dominated by blueshifted absorption from optically 
thick gas flowing outward (see H2.3[l . Instead of the 
high-ionization-state emission lines typically associated 
with symbiotic-star nebulae, this first FU SE observa- 
tion shows a large amount of cool gas (e.g., singly and 



doubly ionized species of C, Fe, and Si). The hot WD 
surface thus appears to have been hidden by an optically 
thick shell of material. During the third stage of the rise 
to optical maximum, Thot slowly began to increase again 
and the U — B color became more blue, the final increase 
in optical flux became dominated by [/-band light, the 
FUV spectra began to move from absorption to emission, 
and the radio flux rose. After the third stage, the X-ray 
spectrum showed evidence for the presence of significant 
quantities of shock-heated gas. 

The uncharacteristically low radio flux density mea- 
sured during the second stage of the outburst (see H2.4[> 
is consistent with the ejection of an optically thick 
shell from the WD. The radio flux is primarily due to 
bremsstrahlung emission from the nebula. A low value 
indicates that the size of the emitting region is unusually 
small, in this case presumably because the flux of ioniz- 
ing photons decreased when the shell first expanded and 
Thot dropped. At the time of the first radio observation 
(Rl), we estimate T ho t < 89, 000 K. We thus expect the 
flux of ionizing photons from the hot-component surface 
to drop more than 80% from the pre-outburst value to 
nioni z {Rl) = 0.16 n i0 niz,q, where n lon i Z {Rl) is the flux of 
ionizing photons from the hot-component photosphere at 



the time of Rl , and : 



Fj wmz,q 



is the flux of ionizing photons 



in quiescence, when Thot ~ 150,000 K. 

Given a typical quiescent 5 GHz radio flux density at 
the orbit al phase corre sponding to Rl of 5 5 GHz.q = 0.9 ± 
0.3 mJy ()Kennvlll995j1 . and a measured 5 GHz radio flux 
density of S5 GHz = 0.42±0.06 mJy, we can infer the 
approximate radius of the hot-component photosphere or 
ejected shell: 



Rhot(Rl) ~Rhot, 



~>5GHz 



(Rl) 



S5GHz,q 

0.17 ± 0.06 Rq, 



1/2 



Tlioniz,q \ 
)j ioniz (^l) / 



1/2 



(4) 



where we have taken the quiescent- state hot-componen t 
radius to be Rhot,q ~ 0.1 R & (see iMurset et all 11991^ . 
We have assumed that changes in the size of the ionized 
portion of the nebula occur in the outer regions that 
are closer to optically thin than optically thick. Thus, 
the change in radio flux density is approximately equal 
to the change in total flux of ionizing photons. Since 
the geometry and ionization structure of the nebula are 
poorly known, the resulting uncertainty in our Rhot es- 
timates from the radio data is probably even larger than 
the quoted errors. We expect the change in flux of ion- 
izing photons to impact the radio emission almost im- 
mediately, since the recombination time in a nebula with 
electron density on the order of 10 9 -10 10 cm -3 is minutes 
to hours l)Fernandez-Gastro et alJlT995[) . 

A similar analysis of the radio flux at the time of 
our second radio observation (R2) gives Rhot(R^) = 
0.23 ± 0.09 R q at JD 2451865, roughly consistent with 
Rhot derived from a FUSE observation taken at the 
same time (see >!3.3l below1. From analysis of additional 
FUSE observations taken between the second and third 
radio observations, we find that Rhot reached a maximum 
size of greater than 0.35 Rq. The third radio observation 
indicates that the hot-component photosphere had con- 
tracted, to Rhot(R3) = 0.27 ± 0.1 R Q , by JD 2451921. 
Thus, the expansion and dissipation of the photosphere 
or optically thick shell is evident in the radio as well as 
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Fig. 10. — [/-band (blue triangles) and V-band (black dots and crosses) light curves from the early part of the 2000-2002 outburst of 
Z And, with the 1997 V-band light curve (red squares) shifted in time and overplotted. The optical rise initially follows the same course 
for the two events. In 2000, however, the first-stage rise proceeds beyond the maximum level reached in the 1997 event. The outburst 
evolution subsequently takes a different path. The blue triangles and black crosses are data from KAIT, and the black dots and red squares 
are data from the AAVSO. 



FUV data. 

Our first radio data point was obtained during the op- 
tical plateau after the second optical rise, over one month 
into the outburst. Therefore, the only information about 
the evolution of Rhot very early in the outburst is that 
indirectly inferred from the decline in That between opti- 
cal spectroscopic observations 01 and 03. A radio light 
curve (plus simultaneous optical spectroscopic monitor- 
ing) starting earlier in a classical symbiotic-star outburst 
would provide more complete information about the evo- 
lution of Rhot during the beginning of the shell-ejection 
phase of this type of event. Since the FUV fluxes can 
provide an independent measure of Rhot , a simultaneous 
FUV light curve early in an outburst would provide a 
test of the shell-ejection hypothesis. 

The non-detection of the magnetic hot spots on the 
surface of the WD at the beginning of the third stage 
of the outburst is also consistent with the ejection of an 
optically thick shell from the WD. Because the WD in 
Z And is magnetic, the accretion flow close to the WD 
is channeled along magnetic field lines. Hot spots at the 
base of these column s produce an optical oscilla tion at 
the WD spin period l|Sokoloski fc Bildstenlll999l) . How- 
ever, on JD 2451872 (just before the optical maximum), 
the oscillation amplitude upper limit was well below the 



amplitude measured during the 1997 outburst, suggest- 
ing that the white-dwarf surface was either hidden or 
the magnetic accretion flow disrupted during the second 
stage of the rise to optical maximum in late 2000. 

Our X-ray observations are also consistent with the 
shell-ejection hypothesis. The hard (>2 keV) tail seen 
in the second X-ray observation (X2) could be the result 
of shock-heated gas, which in symbiotic s tars has been 
sugge sted to be due to colliding winds (e.g.. lMiirset et all 
1997). A new burst of material blown from the WD 
would tend to enhance such emission. Furthermore, as 
discussed in H2.6I at least one of the absorption edges 
wc detected can be ascribed to an ionization state of N 
that is most likely produced in shock- heated gas. The 
difference between the 1997 and 2000 peak X-ray fluxes 
provides additional support for shell ejection during the 
2000-2002 event. In 1997, no shell was ejected, and so 
the accreting WD was clearly visible. Accordingly, the 
X-ray flux reached a significantly higher value in 1997 
compared to 2000, when much of the X-ray flux from 
accretion was probably hidden. 

The physical changes that comprise the third stage of 
the rise are, however, reflected most clearly in the FUV 
spectra. During this period, the FUV emission begins to 
leak through the ejected shell, producing the high FUV 
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fluxes seen in 2000 November and December and the U- 
dominated final optical rise (via nebular reprocessing of 
FUV emission into the optical). The evolution of the 
FUV spectrum from an optically thick, absorption-line 
spectrum with blucshifted P-Cygni absorption to opti- 
cally thin emission-line spectra provides the most direct 
evidence for shell ejection and dissipation. Moreover, 
the evolution of the blucshifted P-Cygni absorption re- 
veals that the average velocity of the absorbing material 
along the line of sight decreased with time. Finally, the 
maximum velocities from the first FU SE observation are 
comparable to the 400 km s _1 inferred for the speed of 
the radio jet that w as produced in conjunct ion with the 
2000-2002 outburst ijBrocksopp et alJl2004j) . 

3.5. Optical Rebrightening in 2002 

In mid-2002, the decline to optical quiescence was in- 
terrupted by a ~1 mag rebrightening that lasted for 
about five months. The U — B color during this re- 
brightening was initially quite red (comparable to the 
second stage of the outburst; see bottom panel of Fig.GJ. 
As the rebrightening event faded, the U — B color be- 
came slightly more blue, although never as blue as dur- 
ing the entire previous year-long optical decline. Both 
the FU SE and VLA fluxes were unusually low during 
the optical rebrightening. 

Our data indicate that the optical rise in 2002 could 
have resulted from a decrease in Thot due to a sli g ht ex- 
pansion of the white-dwarf photosphere. iSkopall (2003) 
proposes that the rebrightening was due to the emergence 
of the hot component from eclipse. The unusually low 
values of the FUV and radio fluxes near the maximum of 
the optical rebrightening, and the orbital inclina tion of 
47 ± 12° from polarimetry l|Schmid fc Schildlll997l) . both 
present problems for the eclipse interpretation. More- 
over, the fact that the optical flux is higher after the dip 
in the light curve despite the system being at an orbital 
phase where it is normally low, and that the minimum of 
the possible eclipse occurs before orbital phase zero (see 
F ig. Et, also tend to favor a non-eclipse explanation. 

Wcbbinkl |)197(tD suggested that a similar rebrighten- 
ing in T CrB was due to the emptying of the accretion 
disk after a pulse of material ejected from the red giant 
formed a ring and then disk around the accreting star. 
For Z And, however, the disappearance of the Ne VI 
lines in FUSE spectrum F8, the reddening of the U — B 
color, and the decrease in the radio flux density all indi- 
cate that That decreased (from around 150,000 K to ap- 
proximately 120,000 K) when the optical flux increased. 
To produce this decrease in Thot, we infer that the WD 
radius probably expanded slightly, despite the fact that 
our FUV data do not allow a reliable estimate of Rhot at 
this time (the FUV flux ratio .Fio59/.Piio3 indicates that 
our FUV continuum fluxes in observations F7-F9 must 
have been contaminated by emission and/or absorption, 
so a WD photosphere model could not be fit to derive 
WD radii). iTomov. Tomova. fc Taranoval l|2004fl report a 
similar conclusion regarding the optical activity in 2002. 

4. DISCUSSION 

In the 2000-2002 eruption of Z And, we see evidence 
for phenomena that usually occur in two distinct types 
of CV outbursts (dwarf novae and classical novae). As 
discussed in H3.3I there are multiple similarities between 



the first stage of the 2000 outburst and the 1997 event, 
which was due to an accretion-disk instability (see H3.2JI . 
The 2000-2002 outburst thus appears to have started 
as an accretion-disk instability. After a few weeks to 
months, however, the outburst energetics became dom- 
inated by an increase in nuclear shell burning on the 
surface of the WD. In addition, during the rise to op- 
tical maximum, an optically thick shell of material was 
ejected, causing 7} lot to drop and P-Cygni profiles to ap- 
pear in the FUSE spectra. We suggest that the 2000- 
2002 outburst of Z And was triggered by the influx of 
hydrogen-rich material from a dwarf-nova-like disk in- 
stability, and then powered by the subsequent increase 
in nuclear shell burning on the WD, in a milder version 
of the phenomenon that powers classical novae. The out- 
burst in Z And therefore combines elements of dwarf no- 
vae and classical novae. We refer to this new type of 
eruption as a combination nova. Since the latter stages 
of typical classical symbiotic outbursts are similar to the 
2000-2002 outburst of Z And, and the the recurrence 
times for classical symbiotic outbursts are shorter than 
theoreti cal calculations can obtain for standard shell 
flashes (jSion. Acierno. fc Tomczvi3 ir979). it is possible 
that many classical symbiotic outbursts are combination 
novae. 

The initial correspondence between the 2000 and 1997 
outbursts suggests that the physical trigger mechanism 
was the same for both events. A disk-instability trig- 
ger mechanism can explain the recurrence times (from 
a few years to more than a decade) of the outbursts of 
Z And and other classical symbiotic stars. Furthermore, 
since symbiotics are wide binaries, their accretion disks 
are probably large. For typical accretion rates of ~ 10 ~ 9 
to ~ 10 _8 Mq yr _1 , the disk temperature ranges from 
~ 40, 000 K close to the white dwarf surface to < 1- 
2 x 10 3 K at ~ 3-R©. These conditions encompass the 
range neede d for an unstable disk (e.g., lCannizzdll993l 
I0sakiin"996l and references therein) . Based on a correla- 
tion between the presence of classical symbiotic outbursts 
and either observational indicators that the red giant is 
close to filling its Roche lobe (and therefore that a disk is 
more likely to be present around the WD) or photomet- 
ric oscillatio ns similar to the superhump s in the SU UMa 
class of CVs. lMikolaiewska et al]l)2002ft suggest that the 
presence of an accretion disk drives the eruptions in the 
symbiotic binaries Z And, YY Her, CI Cyg, BF Cyg, and 
AX Per. Our extensive observations of Z And show that 
this assertion is true for this object, and explain how 
the presence of an accretion disk can produce classical 
symbiotic outbursts. 

For the disk-instability trigger mechanism to produce 
an observable increase in the rate of nuclear shell burn- 
ing on the WD, the amount of material added, AMdi, 
must be a significant fraction of the mass of the WD en- 
velope, Mpnv Extrapolating s l ightly from the curves in 
Fig. 8 of iTownslev fc Bildstenl l)2004l) . we expect the en- 
velope mass on a 0.65 M Q WD accreting at a rate of a 
few times 10~ 8 M & yr" 1 to be around 2 x 10~ 5 M Q . For 
the luminous WD in Z And, the envelope mass should 
be even smaller since less material is required to initiate 
a thermonuclear runaway or weak shell flash on a lu- 
minous WD fe.g.. ISion fc Sta.rrfieldllT98fiT) . Observation- 
ally, Fig. 1101 shows that the first stage of the 2000 event 
rose higher than the maximum reached in 1997, appar- 
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ently crossing some threshold in AMr>i+M env needed to 
produce a thermonuclear response. For the 1997 event, 
which we assume was entirely accretion powered, if wc 
take a bolometric correction BC = —2 mag, we find 
AMdi is approximately equal to a few times 1O _7 M0. 
The amount of material accreted in 2000 was probably 
similar or somewhat larger. Thus AMoi / M env ps 0.01- 
0.1 (for AM DI ps 10~ 7 -10~ 6 M©) appears to be sufficient 
to produce enough of an increase in pressure, and hence 
temperature, to trigger a nuclear response in the burning 
shell of Z And. The relatively small amount of accreted 
material that appears to be able to trigger a shell flash, 
and the correspondingly short amount of time between 
flashes, cannot be explained by current models of nu- 
clear burning on the surface of a WD. Episodic accretion 
through an unstable disk could be the additional factor 
that needs to be taken into account. If we take an av- 
era ge accretion rate into the WD disk of 5 x 10 -8 M 
(see[lijima 2002, who finds typical mass transfer rates in 
symbiotic stars of 10~ 7 to 10~ 8 Af Q yr^ 1 ), then it would 
be possible to both fuel the quasi-steady burning and 
accumulate enough material in the disk to trigger a com- 
bination nova roughly every ten years. 

After the addition of fresh fuel, we expect the nuclear 
burning rate to increase on the thermal time scale at the 
base of the envelope, 



CpTAM C P T AM 



th 



L 



hoi 



M 



(5) 



where Cp is the specific heat at constant pressure, T is 
the temperature at the base of the WD envelope, AM is 
the mass of the envelope, and E nuc = 5 x 10 18 erg g _1 
is the amount of energy released by burning one gram 
of hydrogen-rich material. For steady CNO burning, the 
temperature at the base of the envelope is ~ 3 x 10 7 K. 
However, since the CNO burning rate, e nuc , is a very 
sensitive function of T (e nuc tx T 20 ), and this burning 
rate increased by a factor of ten during the 2000-2002 
outburst of Z And, the observed burning enhancement 
can be produced by about a 10% change in the tem- 
perature at the base of the envelope, or a change of less 
than 0.4 x 10 7 K. Thus, the thermal time scale associated 
with the observed degree of change in nuclear burning in 
Z And is 



tth ~ 1 month 



T 



AM 



4 x 10 6 K J V 2 x 10" 5 M, 



M 



5 x 1O- 8 M yr 



T -(6) 



Here, we have taken an envelope mass of 2 x 10 -5 Mq, 
and an average accretion rate of 5 x 10~ 8 Mq yr _1 . 

The second stage of the 2000-2002 outburst occurred 
about 2 weeks after the beginning of the outburst. This 
response time is quite close to the expected nuclear re- 
sponse time given in Eq. |SJ which could in fact be even 
shorter if the total amount of fresh fuel is a signifi- 
cant fraction of the WD envelope and is dumped onto 
the WD quickly (i.e., within a few days to a week), so 
that the compression becomes adiabatic. To test this 
disk-instability-trigger idea for future symbiotic-star out- 
bursts, it will be important to obtain optical spectra dur- 
ing the first few days of the initial rise. 



Both symbiotic stars and supersoft X-ray sources 
probably contain accreting WDs with nuclear shell 
burning (van den Heuvel et al .111 9921 iMurset et alJl!991t 
iSokoloski et alJl2001^ . However, based upon the lower 
hot-component luminosity estimates for symbiotics, the 
burning may occur at a slower rate in symbiotics (and 
in a few, it is completely absent). The cause of optical 
and X-ray brightness changes appears to be another im- 
portant difference between these two classes of objects, 
both of which are considered to be possible progenitors 
of Type la supernovae. 

A promising model for the long-term variability of 
at least some supersoft X-ray sources is the simple 
photospheric-expansion model (in which the WD bolo- 
metric luminosity remains constant), since X-ray fluxes 
have been seen to drop wh en the optical fluxes ris e in sev- 
eral supersoft sources (e.g.- lSouthwell et al.ll996|) . We do 
not observe this inverse relationship between optical and 
X-ray flux during either the 1997 or the 2000-2002 out- 
bursts of Z And. As discussed in q3.1l and i!3.3l the rise in 
Thot in 1997 and the increase in L^ot by an order of mag- 
nitude during the 2000-2002 outburst are also at odds 
with the simple photospheric-expansion model. Finally, 
given the low (compared to supersoft X-ray sources) hot- 
component effective temperature for Z And, blackbody 
emission from the WD surface probably plays very little 
role in producing the X-ray emission from this source. 
Instead, the X-ray observations of Z And are consistent 
with the combination nova picture. 

5. SUMMARY 

We have shown that the 2000-2002 outburst of Z And 
was primarily powered by an increase in nuclear burn- 
ing on the WD surface. The smaller, 1997 event, on 
the other hand, is consistent with having been a disk 
instability like those of dwarf novae. Thus more than 
one type of classical symbiotic outburst can occur within 
the same symbiotic system. The key points that allow 
us to gain new insight into classical symbiotic outbursts 
from the 2000-2002 eruption of Z And are as follows: 1) 
high-precision UBV light curves that showed the out- 
burst to be a multistage event; 2) the close similarity 
between the first stage of the multi-stage outburst and 
the 1997 outburst, which we identify as a dwarf-nova-like 
disk instability; and 3) the fact that the hot-component 
luminosity during much of the eruption, as determined 
from FUV fluxes, was too high for the outburst to have 
been solely accretion-driven. 

In the 2000-2002 event, the bolometric luminosity in- 
creased during the outburst to ~ 10 4 L Q , indicating that 
the outburst could not have been due to either a sim- 
ple expansion of the white-dwarf photosphere (e.g., in 
response to a small change in M to above the steady- 
burning upper limit) or to a burst of accretion. Instead, 
the overall energetics support thermonuclear shell burn- 
ing as the main source of power. However, since both 
the short recurrence time for large outbursts in Z And 
and the structure of the 2000-2002 outburst light curve 
suggest that the increase in nuclear burning was not a 
standard recurrent nova, we suggest that the enhanced 
nuclear burning was triggered by an accretion event. We 
favor a dwarf-nova-type disk instability as the trigger. 
Thus, we propose a new type of outburst that combines 
the physics of dwarf novae and classical novae. We call 



20 



Sokoloski et al. 



this type of event a combination nova. 
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TABLE 1 
Optical Photometry 



Date 


U 


B 


V 


(JD) 


(mag) 


(mag) 


(mag] 


2451801.75 


9.69 






2451802.75 


9.66 


10.41 


9.39 


2451803.75 


9.76 


10.39 


9.47 


2451804.75 


9.94 


10.61 


9.56 


2451805.75 


9.86 


10.59 


9.53 


2451806.75 


9.76 


10.54 


9.49 


2451807.75 


9.80 


10.57 


9.52 


2451808.75 




10.59 


9.52 


2451810.75 


9.80 


10.44 


9.47 


2451811.75 


9.75 


10.37 


9.46 



Note. — The complete version of 
this table is in the electronic edition 
of the Journal, and includes photome- 
try from 276 nights. The printed edi- 
tion contains only a sample. The rela- 
tive photometry is accurate to approx- 
imately 1%. Absolute photometry is 
accurate to approximately 10%. 



TABLE 2 

Selected Optical Spectroscopic Observations 



Obs. Num UT Date JD-2450000 Days into Outburst 3 



FAST spectra on 162 dates before the 2000-2002 outburst. 



Ol 


2000 Sep 19 


1806 


17 


02 


2000 Sep 23 


1810 


21 


03 


2000 Oct 04 


1821 


32 


04 


2000 Nov 25 


1873 


84 


05 b 


2000 Dec 03 


1882 


92 


06 


2000 Dec 29 


1907 


117 



Plus FAST spectra on 111 additional dates through the end of 2003. 



Note. — See the version of Table |3 in the electronic edition of the 
Journal for the full listing of FAST optical spectra. 

a The beginning of the 2000-2002 outburst was taken to be JD 2451790. 
b From the William Herschel Telescope, courtesy of R,. Hynes. 
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TABLE 3 
Optical Line Equivalent Widths 



Date 
(JD) 


Exposure 
Time (s) 


H 7 
4340A 


[O III] 

4363A 


He I 

4388A 


N III 

4640A 


He II 

4686A 


H/9 
486lA 


He I 

5876A 


[Fe VII] 
6087A 


TiO 
6200A 


Ha 
6563A 


He I 

6678A 


Raman 
6830A 


TiO 
7175A 


2449608.87614 


5.00 


-36.51 


-8.50 


-3.36 


-7.00 


-79.63 


-44.46 


-7.28 


-1.55 


0.59 


-120.30 


-3.96 


-9.69 


1.02 


2449608.87672 


5.00 


-37.96 


-8.25 


-3.21 


-6.42 


-79.45 


-45.29 


-7.27 


-1.23 


0.57 


-123.00 


-3.49 


-9.64 


1.01 


2449608.87730 


60.00 


-38.71 


-9.86 


-4.65 


-7.53 


-81.31 


-45.35 


-7.78 


-1.38 


0.59 


-79.90 


-3.43 


-9.56 


1.01 


2449693.62117 


5.00 


-52.39 


-3.37 


-3.11 


-8.99 


-86.84 


-90.62 


-7.80 


-1.88 


0.55 


-179.72 


-4.18 


-11.34 


0.86 


2449693.62223 


30.00 


-50.95 


-4.87 


-2.97 


-8.63 


-88.37 


-94.11 


-7.32 


-2.12 


0.56 


-171.12 


-4.15 


-11.46 


0.86 


2449725.57784 


5.00 


-51.76 


-2.92 


-3.71 


-8.01 


-84.37 


-96.06 


-7.88 


-1.96 


0.56 


-190.80 


-5.10 


-10.35 


0.84 


2449725.57821 


30.00 


-54.49 


-3.25 


-2.84 


-6.67 


-80.26 


-100.16 


-7.82 


-1.79 


0.58 


-170.60 


-5.03 


-9.94 


0.85 


2449725.57922 


60.00 


-55.38 


-4.08 


-2.60 


-7.41 


-82.00 


-97.73 


-7.93 


-1.81 


0.57 


-153.39 


-5.25 


-10.04 


0.86 


2449749.56896 


60.00 


-55.03 


-3.93 


-2.90 


-6.53 


-78.05 


-95.47 


-10.63 


-1.79 


0.56 


-54.30 


-5.32 


-9.88 


0.76 


2449749.56999 


60.00 


-57.80 


-3.69 


-2.55 


-7.11 


-73.69 


-95.69 


-10.13 


-1.92 


0.56 


-37.69 


-5.38 


-10.98 


0.80 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. The complete version 
contains equivalent widths from 826 spectra between JD 2449608.8761 and JD 2453002.5699. Equivalent widths arc given in units of A and have errors of ~ 
5%-10%. 



TABLE 4 
FUV Observations 



Obs. 


UT Date 


JD 


Obs 


r 


/ a 


/ a 


Num. 




-2450000 


Length (s 


(959.5 A) 


(1058.7 A) 


(1103.4 A) 


Fl 


2000 Nov 16 


1865.1 


8203 


3.6 ±0.2 


7.3 ±0.1 


9.5 ± 0.2 


F2 


2000 Nov 27 


1875.8 


10847 


8.4 ±0.6 


11.1 ±0.1 


13.2 ±0.4 


F3 


2000 Dec 15 


1894.4 


10444 


8.2 ±0.3 


11.1 ±0.3 


13.5 ±0.03 


Fl 


2001 Jul 20 


2111.1 


10928 


3.8 ±0.2 


5.4 ±0.1 


5.8 ± 0.8 


F5 


2001 Jul 22 


2112.8 


13081 


3.8 ±0.2 


5.6 ±0.1 


5.9 ±0.2 


F6 


2001 Sep 30 


2183.1 


10825 


3.7 ±0.3 


5.6 ±0.2 


5.9 ± 0.2 


F7 


2002 Jul 05 


2460.6 


3731 


<0.2 


0.70 ± 0.05 


0.62 ± 0.05 




2002 Oct 20 


2568.5 


59 




Safety Snap 




F8 


2002 Oct 22 


2570.4 


6736 


<0.2 


0.1 ± 0.07 


<0.2 


F9 


2003 Aug 04 


2856.1 


14175 


0.5 ± 0.1 


0.78 ±0.05 


0.72 ±0.07 


a FUSE flux densities in 


units of 10~ 


13 erg cm~ 2 


s" 1 A -1 . 







TABLE 5 
Radio Observations 



Obs. 


UT Date 


Day a 


Instrument 


r 


r 


r 


r 


Num 






(Config b ) 


(1.4 GHz) 


(5 GHz) 


(8.5 GHz) 


(15 GHz) 


Rl 


2000 Oct 13 


1830 


VLA (D) 


< 0.29 


0.42 (0.06) 




0.87 (0.13) 


K2 


2000 Nov 17 


1865 


VLA (A) 


0.28 (0.05) 


0.80 (0.03) 




2.45 (0.09) 


R3 


2001 Jan 12 


1921 


VLA (A) 


< 0.29 


1.21 (0.10) 




2.94 (0.50) 


R4 


2001 Jan 28 


1937 


MERLIN 




1.32 (0.37) 






R5 


2001 Feb 28 


1968 


MERLIN 




0.70 (0.21) 






R6 


2001 Mar 15 


1983 


VLA (B) 


< 0.17 


1.09 (0.03) 




3.50 (0.12) 


K7 


2001 May 08 


2037 


VLA (B) 


0.69 (0.02) 


1.02 (0.02) 




2.53 (0.09) 


R8 


2001 Aug 10 


2131 


VLA (C) 


0.28 (0.08) 


1.06 (0.04) 




2.14 (0.14) 


R9 


2001 Sep 18 


2170 


MERLIN 




1.20 (0.26) d 






R10 


2001 Oct 23 


2205 


MERLIN 




0.86 (0.20) 






Rll 


2002 Apr 30 


2394 


MERLIN 




0.78 (0.20) 






R12 


2002 May 06 


2400 


MERLIN 




0.64 (0.20) 






R13 


2002 Dec 09 


2618 


VLA (C) 




0.50 (0.10) 


0.84 (0.07) 




R14 


2003 Jan 03 


2643 


VLA(C) 




0.42 (0.09) 


0.84 (0.08) 




R15 


2003 Jul 24 


2845 


VLA (A) 




0.65 (0.06) d 




1.55 (0.16) 



a Day = JD-2450000. 

^Config — VLA Configuration. The configuration of the VLA antennas is changed over time, between 
the most extended (A) configuration, with a resolution of about 0.35 arcsec at 4.86 GHz, and the most 
compa ct (D) configuration, with a resolution some 35 times worse. See, e.g., ITavlor. Ulvestad. Perlevl 
I2004D for details. 

c Flux densities in mjy. 

^Spatially resolved. Flux densities quoted for these observations are the total flux densities. 



Sokoloski et al. 



TABLE 6 
X-Ray Observations 



Obs. Num/ 


UT 


Date 


JD 


Orbital 


Obs. Length b 


Source 


Inst 






-2450000 


Phase a 


(ksec) 


Counts 


XI/ Chandra (ACIS-S/HETGS) 


2000 


Nov 13 


1862 


0.12 


19 


109 


X2/XMM - Newton (EPIC-pn) 


2001 


Jan 28 


1938 


0.21 


20/16 


>1500 


X3/XMM - Newton (EPIC-pn) 


2001 


Jun 11 


2072 


0.39 


15/4 


Ril40 



a Orbital phase from ephemcris of Mikolaicwska & Kcnvon ( 1996). 

^Second quantity for the XMM — Newton observations is the remaining observation length after 
removal of high-background data. 



TABLE 7 

Best-Fitting Parameters for the X-ray Spectra 



UT 
Date 


N H 
(10 21 cm" 2 ) 


T kT 

(keV) 


(keV) 


E e dge2 

(keV) 


Xl/dof 


f a 

J ab 


f b 
J unab 


XI (2000 Nov) 
X2 (2001 Jan) 
X3 (2001 June) 


1 4+2-5 

1 S+ 1 - 7 

2 4+ 4 - 8 


21+ 04 

u -' a -2.01 u - ±1 -0.01 

12+ ' 07 
u - ±z -0.07 


64+ 02 


1 02+ - 06 
i ' uz -0.05 

96+ 03 
u - aD -0.03 


0.53/5 
1.11/56 
0.35/11 


2.5 ±0.2 
1.3 ±0.04 
0.2 ±0.05 


5.1 ±0.5 
3.1 ±0.1 
1.4 ±0.3 



Note. — All fit parameters are for blackbody models plus absorption edges, with an additional powerlaw 
component in observation X2 (T is the photon index). The quoted uncertainties are 90% confidence. 



a 0.3-7 keV absorbed flux (10~ 13 erg cm -2 s i). 
b 0.3— 7 keV unabsorbed flux (10 — 13 erg cm -2 s^ 1 ). 



TABLE 8 

Hot-Component Luminosities from FUV Fluxes 



Obs. 


UT Date 


JD 


That 


rp a 


77 a 
-TA1103 


Rhot 


Lhot 


Num. 




-2450000 


(10 3 K) 






(Re) 


(Lq) 


Fl 


2000 Nov 16 


1865.1 


92±20 


250 


240 


0.32 ±0.05 


6600±1700 


F2 


2000 Nov 27 


1875.8 


94±20 


380 


330 


0.36 ±0.06 


9400±2500 


F3 


2000 Dec 15 


1894.4 


95±20 


380 


340 


0.36 ±0.06 


9800±2500 


Fl 


2001 Jul 20 


2111.1 


120±20 


180 


150 


0.18 ±0.03 


6100±1500 


F5 


2001 Jul 22 


2112.8 


120±20 


190 


150 


0.18 ±0.03 


6100±1500 


F6 


2001 Sep 30 


2183.1 


120±20 


190 


150 


0.17 ±0.03 


5600±1400 


F7 b 


2002 Jul 5 


2460.6 


150±30 


24 


16 


0.06 ±0.01 


1500±400 


F8 b 


2002 Oct 22 


2570.4 


120±25 


3 


<5 


0.03 ±0.01 


200±50 


F9 b 


2003 Aug 04 


2856.1 


160±35 


26 


18 


0.06 ±0.01 


2000±500 



"■FUSE extinction-corrected flux densities at 1058.7 A and 1103.4 A in units of 10 -13 
erg cm~ 2 s —1 A -1 . Fluxes are corrected for interstellar extinction using the extinction 
curve of Cardclli, Clavton, & Matins (1989) as implemented in the IDL routine 'unred_ccm' 
(see http:/ /archive. stsci.edu/pub/iue/software/iuedac/procedures/unred_ccm.prol, with E(B — 
V) = U.y 1 / mag and R = A(V)/E(B — V ) = 3.1. 

b For observations F7— F9, the flux ratios ^1059/^1103 are unphysical for a WD photosphere. 
They have therefore likely been contaminated by non-photospheric absorption and/or emission. 
The formal values of Rhot an d L^ot for these observations should thus be used with extreme 
caution. 



